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TO ALL WHOM IT MAY CONCERN: 



Be it known that we, Junhyong Kim, a citizen of the Republic of Korea residing 
at 94 Swarthmore St., Hamden, CT 06517, and Shan Jiang, a citizen of the Peoples 
Republic of China residing at 44 Hillspoint Rd., Trumbull, CT 06611, have invented 
new and useful improvements in 

VISUALIZATION AND MANIPULATION OF BIOMOLECULAR 
RELATIONSHIPS USING GRAPH OPERATORS 



for which the following is a specification. 



VISUALIZATION AND MANIPULATION OF BIOMOLECULAR 
RELATIONSHIPS USING GRAPH OPERATORS 
CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims benefit of U.S. Provisional Application No. 60/221,707, 
filed July 31, 2000. Application Serial No. 60/221,707, filed July 31, 2000, is hereby 
incorporated herein by reference. 

FIELD OF THE INVENTION 

The disclosed invention is generally in the field of analysis of biological 
relationships, and more specifically in the field of computational algorithms for 
representing and analyzing large and heterogeneous molecular biological data. 

BACKGROUND OF THE INVENTION 

Genomics technology has become one of the main driving forces behind 
biomedical research. Information from genomics technology is increasing at an 
exponential pace. Simultaneously, the development of new technologies such as DNA 
microarrays, those of functional genomics, and automatic text retrieval, is greatly 
enriching the kinds of information available. The integration of gene expression data, 
sequence data, and genome annotation would greatly facilitate the utilization of 
genomics information by academic and commercial biotechnology enterprises. 
Accordingly, the synthesis and integration of these disparate sources of genomics data 
into a biologically meaningful information is an immediate and fundamental need. 

Some sources of genomics information such as metabolic pathways traditionally 
are represented in graph form, where nodes or vertices represent genes, and edges or 
arrows represent some biological action between the genes. For example, the Enzyme 
Classification system is a hierarchical graph of enzymes related to each other by 
biochemical action. Other types of information, such as gene function classification, 
have implied graph relationships also. 

However, new genomics technologies such as DNA microarrays are generating 
complex data with no canonical methods of analysis. Complexity in data derived from 
this technology results from both the extreme scale of the data (thousands of 
dimensions) and the uncertainty of the biological implications of measurements such as 
global gene expression levels. Thus a multi-pronged approach to data analysis using 
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various statistical techniques and databases is required in order to achieve a synthesis of 
information. 

The analysis of microarray gene expression data requires the clustering of genes 
into groups of comparable expression profiles across experiments, or the clustering of 
experiments into groups of similar expression patterns across genes. Hierarchical 
clustering (Eisen et al., Cluster analysis and display of genome-wide expression 
patterns. Proc Natl Acad Sci U S A 95: 14863-8 (1998)) and self-organizing maps 
(SOM) (Tamayo et al. (1999) Interpreting patterns of gene expression with self- 
organizing maps: methods and application to hematopoietic differentiation. Proc. Natl. 
Acad. Sci. USA, 96:2907-2912) currently are the algorithms used most commonly for 
expression data clustering, and are implemented in a number of shareware and 
commercial software products. The most salient disadvantage of hierarchical clustering 
is that each individual gene occupies a unique position in the hierarchical tree, and 
cannot be assigned to more than one group. The SOM algorithm requires an arbitrary 
predetermination of the number of clusters to be formed, and thus may yield clusters of 
suboptimal quality. 

In order to overcome the disadvantages of conventional algorithms, several new 
algorithms based on graph theoretic tools have been proposed recently. Ben-Dor et al. 
(1999) Clustering gene expression patterns. J. Comput. Biol., 6(3/4): 281-297, describe 
a clustering algorithm using graph theoretic framework in combination with a 
probabilistic model. They devised an algorithm to generate a clique graph from the 
similarity matrix derived from gene expression data. Input data are represented in a 
disconnected undirected graph in which each gene corresponds to a vertex. A clique 
graph, defined as a disjoint union of complete graphs, represents a possible clustering 
of vertices. This algorithm produces nonhierarchical clusters, the number of which is 
determined by the probabilistic algorithm. 

Another algorithm for expression data clustering was proposed by Sharan and 
Shamir, (2000) CLICK: A clustering algorithm with applications to gene expression 
analysis. ISMB 2000, 307-316, using the graph representation and a statistical model. 
As in the algorithm elaborated by Ben-Dor et al (1999), data elements are represented 
by vertices of a graph. The computation starts from a complete graph, and generates 



multiple subgraphs/clusters by recursively cutting each edge whose weight falls into the 
statistically non-connected category. 

The third algorithm based on graph theory for analyzing expression data, 
biclustering, was developed by Cheng and Church, (2000) Biclustering of expression 
data. IS MB 2000, 93-103. In this algorithm, genes and experiments are represented as 
vertices of a bipartite graph, and are clustered simultaneously. The mean square 
residue score of the data matrix for each cluster is used as a measurement of the 
coherence of gene expression across experiments. The algorithm is designed to find a 
maximum complete bipartite sub-graph with the lowest mean square residue score. 
The result of this computation is a set of gene-experiment clusters in which the 
expression of the genes is coherent across the experiments. Thus, the biclustering 
algorithm creates multiple overlapping clusters that better represent genes that 
participate in multiple pathways. 

Although the algorithms summarized above provide solutions for primary data 
analysis, they do not address the need for comparison, integration, and data mining of 
multiple disparate genomic data sets. To address this need, some data integration 
efforts such as KEGG (Kyoto Encyclopedia of Genes and Genomes) (Kanehisa and 
Goto, (2000) KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids 
Research, 28:27-30; Ogata et al. (1998) Analysis of binary relations and hierarchies of 
enzymes in the metabolic pathways. Biosystems, 47: 119-128; Kanehisa et al. (2000) 
Functional enzyme clusters. Nucleic Acids Research, 28:27-30) and DIP (The Database 
of Interacting Proteins) (Marcotte et al. (1999) A combined algorithm for genome wide 
prediction of protein function. Nature, 402: 83-86; Xenarios et al. (2000) DIP: the 
database of interacting proteins. Nucleic Acids Research, 28:289-91) databases have 
endeavored to integrate into pathways gene relationships previously expressed in binary 
form. However, the computations in these systems were carried out at the database 
level by querying a database for all potential consecutive binary gene pairs, and 
subsequently, integrating them into pathways. Computations carried out within the 
database framework are limited to some relatively simple analyses such as the 
generation of pathways, and coloring genes in the pathway. More complex analyses 
such as comparing disparate data sets, exploring gene network structures, and inferring 



pathways and gene functions, are either beyond the capacity of these systems or 
computationally too expensive to perform. 

BRIEF SUMMARY OF THE INVENTION 

Disclosed is a method for universal representation and integration of 
heterogeneous molecular biological relationships using graph theoretic tools. The 
disclosed invention relates to an electronic system, computer-implemented method, and 
program product in which graphs are stored, manipulated and/or graphically output on 
a display or other output device. Biological molecules are represented as vertices in the 
disclosed graphs. Edges that connect vertices in the graph represent the presence of 
relationships between the molecules. The edge weight of the edges contains 
quantitative or qualitative descriptions of the relationship. Thus, molecular biological 
data of different sources and natures can be represented under a single unified structure 
that provides the foundation for integration of disparate molecular biological data. 
Figure 1 exemplifies the basic components of the disclosed molecular relational graphs. 
Moreover, a complete suite of abstract operations and associated rules are defined for 
the graph such that any specific computation of the disclosed method can be achieved 
by compounding operations according to the rules. Thus operations and rules defined 
for the graph confer powerful tools for assimilating disparate molecular biological data. 

The disclosed method relates to the application of graph theoretical data 
representation coupled with graph operators to biomolecule data analysis. This analysis 
framework is referred to herein as the "molecular relational graphing" (MRG) data 
model or as the "gene-graph operator" (GGO) data model. Using the MRG model, 
analysis techniques for synthesis of disparate sources of knowledge such as those of 
microarray gene expression, protein-protein interaction, and gene function can be 
developed. In some embodiments, the disclosed method relates to the application of 
graph theoretical data representation coupled with graph operators to genomic data 
analysis. 

It is an object of the present invention to provide a system for analyzing and 
graphically visualizing genomic data. 

It is another object of the present invention to provide a comprehensive model 
to organize and store gene relationship information as graphs. 



It is another object of the present invention to provide algorithms to analyze and 
compare molecular relational graphs. 

It is another object of the present invention to provide a software program to 
implement a molecular relational graphing data model. 

It is another object of the present invention to provide a software program to 
visualize the molecular relational graph data. 

It is another object of the present invention to provide a large database for the 
storage and organization of molecular relational graphing data. 

It is another object of the present invention to provide an integrative user 
operation environment based on a graphical flowchart metaphor. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a diagram showing an example of the basic structure of the disclosed 

graphs. 

Figure 2 shows a gene-graph (or molecular relational graph) of protein-protein 
interactions in yeast. Data were generated by yeast two-hybrid assay (Uetz et al., 
2000). Each gene is represented as an oval and the interactions between two genes is 
represented by the line connecting the two ovals. This graph encompassed 1,004 genes 
and 957 interactions. Approximately 500 genes form the largest interconnected 
structure. The rest form a number of smaller structures. 

Figure 3 shows a gene-graph (or molecular relational graph) of gene ontology 
functional relationships for a selected set of yeast genes. Thirty-one genes are included 
in this graph. Their participation in multiple functional processes makes the 
intersecting pathways form a dense network. 

Figure 4 shows a gene-graph (or molecular relational graph) of expression 
analysis data. Data were from a correlation analysis of microarray hybridization 
experiments reported by Spellman et al. (1998). Edges in the graph represent the 
correlation between two genes in gene expression profile. This graph is derived by 
edge-thresholding at 0.4. This graph is generated from correlation analysis of yeast 
gene expression profile during cell cycle. 

Figures 5A, 5B, 5C, 5D, and 5E show a gene-graph analysis (or molecular 
relational graphing analysis) of expression data from microarrays hybridizations assay. 
Figure 5A shows the gene-relationship structure derived by applying the AND operator 



between the Gene Ontology (GO) annotation graph and the gene expression graph, 
wherein both graphs have the same graph structure. Two structures are labeled as *1 
and *2, respectively. Figure 5B shows the expression gene-graph threshold at 0.1. 
Both structure *1 and *2 are present, some relationships are missing in structure *1 due 
to the high-stringency thresholding. One novel structure (V) cannot be derived from 
naive GO annotation grouping. However, it is supported by the sophisticated grouping 
as shown in Figure 5E. Figure 5C shows an expression gene-graph thresholded at 0.2. 
Both structure *1 and *2 are completely preserved, and the novel structure V is 
expanded by the addition of one gene and two new relationships. Figure 5D shows an 
expression gene-graph thresholded at 0.3. Structure *1 is completely preserved while 
*2 is expanded into a larger one with additional genes and relationships. Structure V is 
expanded also and a fourth structure appears in the graph. Figure 5e shows the relative 
positions of two GO id numbers GO:0007330 and GO:0007328 in GO annotation tree. 
This GO genealogy clearly indicates the legitimacy of the relationship that forms the 
structure V. 

Figure 6 is a diagram of an overview of an example of the design of a data 
mining system using the disclosed method. 

Figure 7 is a diagram of an example of the design of a data mining service 

client. 

Figure 8 is a diagram of an example of the design of a data mining service 

broker. 

Figure 9 is a diagram of an example of the design of a graph computation 
manager. 

Figure 10 is a diagram of an example of the design of a graph computation 

engine. 

Figure 11 is a diagram of an example of the design of a graph visualization 

engine. 

Figure 12 is a diagram of an example of the design of a graph computational 

library. 

Figure 13 is a diagram of an example of the design of a data interface. 
Figure 14 is a diagram of an example of a general purpose computer 
implementing an example of the disclosed method and composition. 
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Figure 15 shows a Unified Modeling Language diagram of GGO (or MRG) 

objects. 

DETAILED DESCRIPTION OF THE INVENTION 

Disclosed is a method for universal representation and integration of 
heterogeneous molecular biological relationships using graph theoretic tools. In the 
method, biological molecules can be represented as vertices in the graph. Edges that 
connect vertices in the graph can represent relationships between molecules. Edge 
weight can contain quantitative or qualitative descriptions of the relationship. In this 
way, molecular biological data of different sources and natures can be represented 
under a single unified structure that provides the foundation for integration of disparate 
molecular biological data. Moreover, a complete suite of abstract operations and 
associated rules can be defined for, and applied to, the graph such that any specific 
computation of the disclosed method can be achieved by compounding operations 
according to defined and devised rules. Thus, operations and rules defined for the 
graph confer powerful tools for assimilating disparate molecular biological data. 

The disclosed method is referred to herein as molecular relational graphing 
(MRG) and involves generation and manipulation of graphs, referred to herein as 
molecular relational graphs. Alternatively, the method is referred to as gene-graph 
operator (GGO) and the graphs are referred to as gene-graphs. 

The disclosed method can be implemented as computer software. For example, 
a molecular relational graphing software program can be written using any suitable 
programming language, such as the Java™ programming language. A software 
program implementing the disclosed method can have two principal features: (1) 
implementation of molecular relational graphing objects and the ability to store in a 
local and/or remote database, and (2) implementation of operators. Such operators 
manipulate the molecular relational graphs as objects, much as mathematical operators 
manipulate numbers. Like mathematical operators, molecular relational graphing 
operators allow direct manipulation of graphs using graph operations such as addition 
and subtraction. 

Molecular relational graphing is preferably implemented on a programmed 
general purpose computer system. However, the molecular relational graphing can also 
be implemented on a special purpose computer, a programmed microprocessor or 



microcontroller and peripheral integrated circuit elements, an ASIC or other integrated 
circuit, a hardwired electronic or logic circuit such as a discrete element circuit, a 
programmable logic device such as a PLD, PLA, FPGA or PAL, or the like. 

The disclosed molecular relational graphing method provides a comprehensive 
framework to accommodate disparate data sets; the underlying graph theoretic tools 
confer powerful approaches, for example, to analyze network structures, and to infer 
pathways and functions. The method complements existing integrative efforts. Most 
importantly, the integrative and analytical capacity of the disclosed molecular relational 
graphing is far greater than that of any existing algorithm. 

The disclosed method provides a new technique for genomics data analysis, 
including that generated by microarrays. In the disclosed method, heterogeneous 
genomics information can be unified into a common graph-theoretic structure. 
Subsequently, formal graph operators can be defined, allowing the manipulation of 
different information through a syntax of graph structures. The disclosed method allows 
querying of complex information with a dynamic rearrangement and synthesis of 
heterogeneous data. 

The disclosed method offers a universal representation of heterogeneous 
molecular biological data. Biological data of different sources can be captured in a 
single unified structure based on intermolecular relationships. Modification and 
integration of heterogeneous data are achieved by applying single or compounded 
operations on multiple data sets. Thus, unlike previous techniques, the disclosed 
method is not restricted to any particular problem domain and is not limited to a few 
fixed kinds of data integration. As used herein, heterogeneous biological data, 
heterogeneous molecular biological data, or heterogeneous biomolecular data refers to 
data from different types of biological systems (thus embodying different types of 
relationships between biological molecules), different types of measurements (thus 
embodying different types of relationships between biological molecules), different 
types of biological molecules (preferably different types of biological molecules that 
have relationship with each other), or any other combination of disparate biological 
data. As an example, one form of heterogeneous molecular biological data would be 
expression relationships between genes and proteins (two different types of biological 
molecules). Another form of heterogeneous molecular biological data would be the 
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combination of a variety of expression and physiological measurements (that is, 
multiple different relationship nd biological molecules) for a particular type of cell or 
tissue. 

Different types of biological systems include, for example, protein-protein 
interactions; protein-nucleic acid interactions; gene expression regulation; protein 
expression regulation; cellular signal transduction pathways; physiological states; 
disease states; and metabolic pathways. Different types of measurements include, for 
example, the presence of association in time, or space, or logical meaning; physical or 
logical states such as activation and inhibition; real value measurement of spatial 
distance such as physical distances between genes, open reading frames, single 
nucleotide polymorphisms, expressed sequence tags, sequence tag sites, or a 
combination thereof; sequence similarity between genes or proteins; structural 
similarity between proteins; radiation hybrid mapping distances between genes, open 
reading frames, single nucleotide polymorphisms, expressed sequence tags, sequence 
tag sites, or a combination thereof; genetic distances between genes, open reading 
frames, single nucleotide polymorphisms, expressed sequence tags, sequence tag sites, 
or a combination thereof; real value measurement of time or kinetic information such as 
chemical conversion rate; Euclidean and other distance metrics in feature space to 
measure logical relationship; correlation coefficient as a statistical metric to measure 
logical relationship; values of fuzzy set membership function as a metric to measure 
logical relationship; and conditional probability as a measurement of causal 
relationship. 

Different types of biological molecules include, for example, genes, open 
reading frames, expressed sequence tags, single nucleotide polymorphisms, sequence 
tag sites, nucleic acids, DNA, RNA, mRNA, cDNA, proteins, peptides, enzymes, 
metabolites, carbohydrates, exons, introns, cleavage fragments, restriction fragments, 
amino acid modifications, protein domains, DNA or RNA secondary or tertiary 
structures, nucleic acid motifs, protein motifs, and metal ions. 

In the context of the disclosed molecular relational graphs, use of heterogeneous 
molecular biological data is manifested by having at least two of the vertices represent 
different types of biological molecules; having at least two edges represent different 
types of relationships between the biological molecules represented by the vertices 



connected by the edges; having at least one edge represent a plurality of different types 
of relationships between the biological molecules represented by the vertices connected 
by the edge; and/or having at least one vertex represent a plurality of different types of 
biological molecules. 

A graph is a mathematical abstraction of relationships among different entities 
in the real world. The graph represents an entity (such as a gene, protein, or other 
biomolecule) as a vertex, and encapsulates the relationship between two entities as an 
edge that connects the two vertices. The interconnections among a set of vertices, 
designated by a set of edges, form a graph. Many algorithms have been developed that 
allow efficient manipulation of the graph, retrieval of information stored in the graph, 
and computation using graphs as objects. Graph theory and techniques can be applied, 
in the disclosed method, to model and manipulate biomolecules and biological 
relationships organized as a graph. 

The disclosed method relates, in part, to the application of the gene-graph 
operator method to the analysis of genomic relationships. Genomic relationships can 
be encapsulated by a graph model regardless of the context and the technology from 
which the information is derived. In GGO, each gene (or protein or biomolecule) is 
represented as a vertex in the graph, and the relationship between two genes (or 
proteins or biomolecules) is represented as the edge between vertices. The graph model 
can be used to represent various types of genomic relationships (or other biomolecular 
relationships) as defined by the contents of the vertex and the edge. For example, a 
graph can model a gene expression data set if the edge contains the measurement of 
correlation of the expression patterns of two genes. With such a gene-graph model, 
algorithms developed in graph theory enable sophisticated analysis of the gene- 
relationship data. Examples of complex analysis include the elucidation of mechanisms 
of gene regulation, the identification of gene action pathways, and the identification of 
critical genes that link multiple biochemical pathways. 

In some embodiments, the disclosed method can use and manipulate large 
databases, including object-oriented databases, for the storage and organization of 
molecular relational graph data (or gene-graph data), and can implement molecular 
relational graphing models for proteome and genome mapping data. A molecular 
relational graphing database can comprise large data sets from a variety of sources, 
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such as gene expression analysis, proteome analysis, genome mapping, and functional 
genome annotation. Data objects, n-nary operations, and graph functions can be 
implemented as, for example, individual software components, which then can be 
connected to implement a particular set of analysis operations. The software 
components can be graphically represented as iconized tools. Connections between 
components can be established by the user from a graphical interface. 

The manipulations of graphs in the disclosed method may involve single graphs 
(by using unary operators) or multiple graphs (by using binary and n-nary operators), 
and may produce numerical results or new graphs (referred to herein as product 
graphs). These manipulations can be designed such that they can be combined into a 
sequence of steps to produce a particular synthetic meta-analysis. The manipulations 
can also be recursive, with, for example, a result of a manipulation being manipulated 
again (or multiple times) in the same way. The results of the meta-analysis can be 
interpreted in a biological context. In other words, instead of fixing the results of, for 
example, microarray analyses or various genomics information into a static and 
awkward data model, the information can be encapsulated into a common graph 
structure with associated syntactic rules that are defined for manipulating the common 
structure. This encapsulation produces an information model that is dynamic and 
particularly suited to synthesis of disparate information. 

The disclosed method and composition can be understood further by reference 
to the following example system, which describes an example of the use of a gene 
graph operator (which is also referred to as a molecular relational graphing operator) at 
the heart of a data mining and interface system. The gene graph operator (Figure 12) is 
a software embodiment of the disclosed method and provides representations for all 
types molecular relational graphs (gene-graphs). The gene graph operator is used by 
the graph computation executor in the graph computation engine (Figure 10) to 
construct molecular relational graphs and perform operations on molecular relational 
graphs. 

As illustrated in Figure 6, the user can submit a data mining request by 
interfacing with the data mining service client (details in Figure 7). The data mining 
service client includes the user interface and displays results of data mining and graph 
manipulation (Figure 7). The data mining service client then makes a data mining 

11 



request of the data mining service broker (details in Figure 8). The data mining service 
broker decomposes data mining requests and dispatches requests for data to various 
subsystems. The data mining service broker also communicates the results of data 
mining, graph construction, and graph manipulation to the data mining service client. 
5 As illustrated in Figure 6, the data mining service broker makes graph 

computation requests to the graph computation manager (Figure 9). The data mining 
services broker also receives the results of data mining, graph construction, and graph 
manipulation from the graph computation manager (Figure 6). The graph computation 
manager interfaces with databases to receive graph data (Figure 6). The graph 
10 computation manager sends graph computation requests to the graph computation 
engine (Figure 10). The graph computation engine builds graphs from the data 
q received from the graph computation manager and performs operations on graphs. The 

results of the computations are communicated to the graph computation manager 
*** (Figure 6). The graph computation manager also sends graph visualization requests to 

g 15 the graph visualization engine (Figure 11). The graph visualization engine produces 

graphics objects from graph data and communicates the graphics objects to the graph 
computation manager (Figure 6). The graph computation manager sends the graphics 
objects and non-graph data from data mining operations to the data mining service 
broker which in turn communicates the non-graph data and graphics objects to the data 
20 mining service client where the user can access and view the results (Figure 6). 

The disclosed method and composition can be understood further by reference 
to the following example system. As illustrated in Figure 14, the user can load data and 
interact with the system through network interface 110, disk 118 and 114, keyboard 
124, or a combination. The user graph data can be formatted as flat files of ASCII or 
25 binary type; files with fields separated by comma, tab, line break, carriage return, or 
paragraph or other character codes for import into spreadsheets. A preferred format is 
appropriate tables of a relational database. The graph data can be accessed by a graph 
manipulation component such as GGO subsystem 102 (see also Figure 6). The GGO 
subsystem can obtain graph data by request from the data mining service broker 104 
30 (see also Figure 8). The system can display for the user visual representations of graph 
data on monitor 126 or other display device. 
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To adapt graph structures to the analysis of biomolecule relationship data, graph 
theoretical vocabulary can be defined in a biological context. Using this vocabulary, 
biomolecular relationship information, such as information derived from gene 
expression analysis or the Gene Ontology (GO) database, can be represented and 
integrated using the disclosed molecular relational graphing model. 

Accordingly, for purposes of the disclosed method, by "graph" it is meant a 
collection of vertices (nodes) and edges denoted as G = { V, E} where V is the set of 
vertices and E is the set of edges. 

By "vertex" and "vertices" it is meant an encapsulation representing a biological 
molecule such as DNA, RNA, protein, or small compounds. Vertices can be labeled 
with the identities of the biological molecules. If two different graphs share identically- 
labeled vertices (or one or more allowed aliases), it is assumed, unless the context is to 
the contrary, that they are comparable. For example, a vertex in a gene expression 
graph might be labeled "CDC28" and a vertex in a protein-protein interaction graph 
might also be labeled "CDC28". They are assumed to be comparable even though the 
actual molecules in the experiments might not be identical. Vertices can encapsulate all 
the properties of the biological molecules, and therefore, may be multi-labeled. 

By "hyper- vertex" it is meant a set of vertices representing a set of biological 
molecules. Unless the context clearly indicates otherwise, the term "vertex" is used 
herein to refer to both vertices as defined above and hyper- vertices. 

By "edge" is it meant a connection between two vertices. It usually represents a 
relationship between the biological molecules specified by the two vertices. An edge 
can be directed, representing the direction of action, and it can be weighted. An edge 
can be said to be defined by a pair (a, b) where a and b each represent a vertex. 

By "edge weight" it is meant a number or a descriptor assigned to an edge, 
denoting a quantitative degree of relationship or qualitative type of relationship. For 
example, a real-valued edge weight can denote the correlation coefficient between 
expression patterns of two genes; an edge weight with the descriptor "+" can denote 
"activation" of one gene by another. 

By "hyper-edge" it is meant an edge which connects two or more vertices as a 
set denoting a relationship that involves more than pair-wise interactions. A hyper- 
edge may also be weighted. A hyper-edge can be said to be defined by a pair (a, b) 
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where at least one of a and b represents a set of vertices. For a regular hyper-edge, both 
a and b represent a set of vertices. Unless the context clearly indicates otherwise, the 
term "edge" is used herein to refer to both edges as defined above and hyper-edges. 

By "directed edge" it is meant an edge defined as an ordered pair (a, b) where a 
and b are vertices. 

By "undirected edge it is meant an edge defined as an unordered pair (a, b) 
where a and b are vertices. 

* 

By "directed hyper-edge" it is meant a hyper-edge defined as an ordered pair (a, 
b) where a and/or b are sets of vertices. 

By "undirected hyper-edge it is meant a hyper-edge defined as an unordered 
pair (a, b) where a and/or b are sets of vertices. 

In some embodiments, the disclosed software can perform the task of 
integrating data from, for example, microarray gene expression analysis, Gene 
Ontology annotation, and protein-protein interaction analysis into a molecular relational 
graphing data model. The disclosed software can also have functions for pathway 
analysis, critical gene identification, gene-action subsystem identification, and pathway 
comparison. Since the molecular relational graphing model is best illustrated using a 
graphical approach, also disclosed is visualization software for the demonstration of 
data resulting from computation using the disclosed molecular relational graphing data 
model. Such software can be written in any suitable programming language, for 
example, the Java programming language. 

Graph objects, n-nary operators, and graph operators can be implemented as 
individual software components, which are then connected in series using connectors to 
implement the desired set of analysis operations. The software components and 
connectors can be graphically represented as intuitively recognizable glyphs. The user 
of the software can establish connections between components by using the graphical 
interface. Standard analysis techniques can be integrated into the disclosed analysis 
platform by incorporating standard commercial software packages. This will allow the 
system to use many analysis features from other packages, such as clustering analysis, 
for preliminary data processing. The resulting data can be transformed into the 
molecular relational graphing model for high-level analysis. 
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In some embodiments, molecular relational graphing models for proteome and 
genome mapping data will be used. In such embodiments, the molecular relational 
graphing database can contain large data sets from gene expression analysis, proteome 
analysis, genome mapping, and/or functional genome annotation. 
A. Graph Elements 

The disclosed method uses graphs to embody and manipulate relationships 
between biomolecules. Heterogeneous molecular biological relationships can be 
effectively encapsulated in different molecular relational graphs. In a molecular 
relational graph, biological molecules are represented by vertices and information of 
relationships between molecules is stored in edges connecting vertices. 

1. Vertices 

Different types of biological molecules can be represented as different types of 
vertices in molecular relational graphs. Biological molecules that can be represented 
by vertices in molecular relational graphs include but are not limited to: 

genes, open reading frames, expressed sequence tags, single nucleotide 
polymorphisms, sequence tag sites, nucleic acids, DNA, RNA, mRNA, cDNA, 
proteins, peptides, enzymes, metabolites, carbohydrates, exons, introns, cleavage 
fragments, restriction fragments, amino acid modifications, protein domains, DNA or 
RNA secondary or tertiary structures, nucleic acid motifs, protein motifs, and metal 
ions. 

As used herein, "biological molecule" and "biomolecule" refer to any molecule 
or portion of a molecule or multi-molecular assembly or composition, that has a 
biological origin, is related to a molecule or portion of a molecule or multi-molecular 
assembly or composition that has a biological origin. Biomolecules can be completely 
artificial molecules that are related to molecules of biological origin. 

The content of a vertex can include a label and an information table. To 
construct a vertex, a name that uniquely labels a biological molecule can be used as the 
label for the vertex. Properties of the biological molecule can be stored in an 
information table as a part of the content possessed by the vertex such that each row of 
the table contains a property name and a property value. 

Using information retrieved from the Sacchoromyces Genome Database (SGD) 
(Cherry et al., Sacchoromyces Genome Database), the following illustrations provide 
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examples of constructing vertices representing yeast open reading frames (ORFs), 
protein molecules, and genes. 

Illustration 1: Defining vertices representing yeast open reading 
frames (ORFs). 

More than 5,000 genes were identified in yeast genome by either experimental 
or computational methods (Cherry et al. (1997)). Each gene consists of one or more 
exons in its genomic sequence that, when spliced together in order, forms the sequence 
of mRNA for this gene. Part of the mRNA molecule will be translated into proteins. 
The translated portion of the mRNA molecule sequence does not contain any 
translational stop codon. Thus, a continuous fragment of genomic sequence, which 
constitutes a part or whole of translated portion of an mRNA molecule, can be named 
an open reading frame (ORF). 

To construct vertices representing yeast ORFs (Cherry et al. (1997)), a unique 
label for a vertex can be specified, for example, using the name of the ORF such as 
"YCL040W". A vertex can also possess an information table in which properties of the 
represented yeast ORF can be stored. The information table can have two columns: 
<property_name> and <value>. The content of the table can comprise a set of 
(property j%ame, value) pairs that can include, for example: alias, 
chromosomejocation, genomic_sequence_source, description, gene_product, function, 
cellular_component, process, and phenotype. Table 1 shows the content and structure 
of the information table for a vertex representing a yeast ORF, YCL040W. 



Table 1. Information table for a vertex representing yeast ORF YCL040W. 



Property_name 


Value 


Alias 


GLK1 


chromosome ^location 


chromosome_3 


genomic_sequence_source 


SGD_YCL040W 


Description 


Glucose phosphorylation 


genejproduct 


Glucokinase 


Function 


Glucokinase 


Cellular_component 


Cytosol 
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Process 


Glycolysis 


Phenotype 


Null mutant is viable with no 
discernible difference from wild- 
type; hxkl, hxk2, glkl triple null 
mutants are unable to erow on any 
sugar except galactose and fail to 
sporulate. 



Illustration 2: Defining vertices representing yeast proteins. 

To represent yeast protein molecules using vertices, one vertex can represent 
one protein molecule. In this representation, the label of a vertex can be assigned the 
name of the represented protein molecule. An information table can be constructed for 
each vertex. The table can comprise two columns: <property_name> and <value>. A 
list of (property _name, value) pairs can be stored in the table. In the information table 
possessed by different vertices, the same property _name may be associated with 
different values. The list of property _names can include, for example: alias, 
sequence_source, structure, EC_number, description, function, cellular_component, 
process, and phenotype. An information table for a vertex representing yeast protein 
grxl is shown in Table 2. The label of the vertex is GRX1. 



Table 2. Information table for a vertex representing yeast protein grxl. 



Property_name 


Value 


sequence_source 1 


PID_G5328 


sequence_source2 


SwissProt_P25373 


sequence_source3 


PIR_S 19363 


Structure 


Sacch3D_YCL035C 


Description 


Glutaredoxin 


Function 


Glutaredoxin 


cellular_component 


Unknown 


Process 


oxidative stress response 


Phenotype 


Null mutant is viable but sensitive 
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Illustration 3: Defining vertices representing yeast genes. 

A complete representation of yeast genes can consist of information for both the 
genomic sequence and the protein products of the gene. By merging together 
5 information contained in vertices representing the ORFs of a gene and the 

corresponding protein products, a vertex that represents the gene can be constructed. 
To create a vertex representing a yeast gene, given that a vertex (vertices) representing 
the ORF(s) of the gene and a vertex (vertices) representing the protein product(s) of the 
gene are created previously, a series of operations can be performed. For example: 
10 Assign the name of the gene to the label for the vertex. 

Create an information table for the vertex. 

Add (property jname, value) pairs (ORF, ORF_name) to the table. ORF_name 
is the label for a merged-in vertex representing an ORF. There may be several (ORF, 
ORF_name) pairs if the gene encompasses more than one ORF. 

15 Add the second type of (property _name, value) pairs, (protein, protein_name), 

to the table. Protein_name is the name of the merged-in vertex representing a protein 
molecule. There may be several (protein, protein_name) pairs if the gene is translated 
into protein molecules of more than one isoform. 

Add additional (property _name, value) pairs to the table such that each pair 

20 consists of the label of a merged-in vertex and the information table possessed by the 
corresponding vertex. 

As an example, a vertex representing a yeast gene, GRX1, is created from a 
vertex representing an ORF, YCL035C, and a vertex representing a protein molecule, 
grxL Since the gene contains only a single ORF and a single protein product, there is 

25 only one ORF vertex and one protein vertex participating in the construction of the 

vertex representing the gene. The label of the vertex representing the gene is specified 
as GRXL The information table for the vertex is shown in Table 3. 
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Table 3. Information table for a vertex representing yeast protein grxl 
Property .name Value 



ORF1 



Protein 



YCL035C 



grxl 



YCL035C 

Property__name Value 



chromosome_location 

Sequence coordination 

genomic_sequence_source 

Description 

gene__product 

Function 

Process 

Phenotype 



chromosome_3 

61173 to 60841 

SGD_YCL035C 

Glutaredoxin 

Glutaredoxin 

Glutaredoxin 

oxidative stress response 

Null mutant is viable but sensitive to oxidative 
stress, grxl grx2 null mutants are viable but lack 
heat-stable oxidoreductase activity. 



GRX1 



Property_name 

sequence_source 1 

sequence_source2 

sequence_source3 

Structure 

Description 

Function 

cellular_component 



Value 

PID_G5328 

SwissProt_P25373 

PIR_S 19363 

Sacch3D_YCL035C 

Glutaredoxin 

Glutaredoxin 

Unknown 
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Process 



oxidative stress response 



Phenotype Null mutant is viable but sensitive to oxidative 

stress, grxl grx2 null mutants are viable but lack 
heat-stable oxidoreductase activity. 

2. Edges 

Information about relationships between biological molecules can be 
represented by edges of molecular relational graphs. Types of quantitative or 
qualitative measurements of relationships stored in edges can include but are not 
limited to the following: 

boolean values indicating the presence of association in time, or space, or 
logical meaning, descriptors of physical or logical states such as "+" representing 
activation and "-" indicating inhibition, real value measurement of spatial distance such 
as physical distance between two genes on the chromosome, real value measurement of 
time or kinetic information such as chemical conversion rate, Euclidean and other 
distance metrics in feature space to measure logical relationship, correlation coefficient 
as a statistical metric to measure logical relationship, values of fuzzy set membership 
function as a metric to measure logical relationship, conditional probability as a 
measurement of causal relationship, and any combination of these. 

Relationships embodied in the disclosed edges can also include physical 
distances between genes, open reading frames, single nucleotide polymorphisms, 
expressed sequence tags, sequence tag sites, or a combination thereof; genetic distances 
between genes, open reading frames, single nucleotide polymorphisms, expressed 
sequence tags, sequence tag sites, or a combination thereof; protein-protein 
interactions; protein-nucleic acid interactions; gene expression regulation; protein 
expression regulation; cellular signal transduction pathways; sequence similarity 
between genes or proteins; structural similarity between proteins; radiation hybrid 
mapping distances between genes, open reading frames, single nucleotide 
polymorphisms, expressed sequence tags, sequence tag sites, or a combination thereof; 
and metabolic pathways. 

The content of an edge can include, for example: (a) labels of two vertices that 

are connected by the edge; (b) directional labels for the two vertices such as "head" and 

"tail" indicating the direction of the edge if the relationship is directional between the 
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two biological molecules represented by the two vertices; and (c) an edge weight table 
which stores properties of the relationship between the two represented biological 
molecules. The edge weight table of an edge can be organized such that each row of 
the table contains a label for a relationship property and a value for the corresponding 
property. 

In the disclosed graphs, vertices represent involved biological molecules and 
edges represent relationships between molecules. Thus the relationship information 
stored in the edge can include, for example, the identities of participating molecules, 
the nature of the relationship, and the properties of the relationship. The following 
illustrations provide examples of creating different types of edges to encapsulate 
different types of relationship information. As used herein, "relationship" refers to any 
characterization shared with, linking, correlating, identifying, or otherwise describing 
any two or more objects (such as biological molecules). 

Illustration 4: Defining edges representing the relationship of protein- 
protein interaction between yeast protein molecules. 

Whole genome-scale study of protein-protein interactions has been carried out 
for yeast (Uetz et al. (2000)). Out of more than 6,000 proteins, 1,004 yeast proteins 
were reported to participate in 957 physical interactions with other protein molecules in 
yeast two-hybrid assays. In order to study large number of protein-protein interactions 
found in yeast cells, interactions between yeast protein molecules can be represented 
effectively using edges defined in molecular relational graphs. To define an edge 
representing a physical interaction between a pair of yeast proteins, vertices 
representing the two participating protein molecules can be defined first. Once the 
vertices are defined, an edge can be defined by, for example, the following three 
components: 

(1) Labels of input vertices and output vertices representing the involved protein 
molecules. 

(2) A Boolean variable, DIRECTED, representing whether the edge is directed 
(thus respecting the input to output designation) or undirected. Since the protein- 
protein interactions are symmetrical relationships for this example, DIRECTED = 
FALSE. 
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(3) An edge weight table in which (property, value) pairs reflecting the 
properties of relationships are stored. In the simplest form, the table contains a list of 
(property, value) pairs such as: (assay_system, two hybrid), (assay_method, beta gal), 
and (strength, 1200). 

Assay_method indicates that the lac-Z gene is used as a reporter and (3- 
galactosidase activity mediates the reporter gene activation and the experimental read- 
out for the assay system. Thus, in this example, the measurement of the strength of 
interaction is a spectrophotometric measurement of absorption of yeast lysate incubated 
with |3-galactosidase substrate. 

To encapsulate the yeast protein-protein interaction data set published by Uetz 
et al. (2000), 1,004 vertices are created to represent all the involved proteins and 957 
edges are created to connect vertices representing the interacting protein pairs. 

Illustration 5: Defining edges representing metabolic pathways in the cell. 

In the cell, metabolic molecules such as glucose and amino acids are 
transformed by various enzymes into different kinds of molecules continuously. These 
metabolites are either disintegrated into simpler molecules or integrated with other 
molecules or modified to form more complex molecules. These pathways of molecular 
transformation can be encapsulated using vertices and edges. To do so, metabolites can 
be represented by vertices first such that each metabolite is represented by one vertex. 
Properties of a metabolite such as the name of the chemical compound, the database 
source of the molecular structure, and cellular localization of the molecule can be 
stored in the vertex. In the representation of metabolic pathways, an edge can be used 
to encapsulate a set of metabolic reactions catalyzed by a given enzyme. Thus, an edge 
connects a pair of vertex groups, one of which represents a group of reaction substrates 
and the other of which represents a group of reaction products. The definition of an 
edge for metabolic pathways can comprise, for example, the following information: 

(1) A set of labels of input vertices representing reaction substrate molecules; 

(2) A set of labels of output vertices representing reaction product molecules; 

(3) DIRECTED = TRUE; 

(4) An edge weight table can be constructed to contain (property __name, value) 
pairs of a list of properties including, for example: 

(a) Enzyme_name: the name of the enzyme that catalyzed the reaction; 
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(b) K m : the Michaelis-Menton reaction rate coefficient; 

(c) Vmax: maximum reaction rate under Michaelis-Menton model. 
Thus, the edge weight table can encompass information about the identity of the 

enzyme that catalyzes the reactions and the kinetics that describe the behaviors of the 
enzyme and the characteristics of the reaction. 

Illustration 6: Defining edges representing functional relationships between 
genes of an organism. 

Functional relationships between genes are summaries of various relationship 
information about the functional roles played by these genes. One example of these 
functional relationships between two genes is that two genes are co-regulated in 
transcription by the same transcriptional factor. Another example is that protein 
products of two genes are immediate neighboring elements in a cellular signal 
transduction pathway. A third example is that protein products of two genes participate 
in the formation of the same holoenzyme complex. Each edge can encapsulate one 
elementary type of functional relationship. Multiplexed complex functional 
relationship representation can be derived using graph operators as discussed below. 

To define edges representing functional relationships between two yeast genes, 
vertices representing the two genes should be defined first. Given the vertices 
available, an edge can be created to represent each elementary type of functional 
relationships between two genes. An edge can be constructed by defining a list of 
information components including, for example: 

(1) Labels of input and output vertices representing the two yeast genes - 
vertex_labell and vertex_label2. 

(2) Assignment to the variable DIRECTED. For example, for signal 
transduction pathways, DIRECTED = TRUE. 

(3) An edge weight table of properties of the elementary type of functional 
relationship stored as (property _name, value) pairs. For example, suppose a protein 
product of gene 2 is a ligand molecule that engages a receptor that is the protein 
product of gene 1 and the ligand-receptor binding activates the next step of signal 
transduction cascade. To represent this type of functional relationship, an edge weight 
table can be constructed to contain (property _name, value) pairs such as: 

(Relationship_type, signal transduction) 
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(Relationship_measurement, Kd) 
(Kd, ligand_binding_constant), 

where K4 is the binding constant which is the measurement of the kinetics of 
binding process. 
B. Graphs 

The disclosed vertices and edges make up the disclosed molecular relational 
graphs. A graph can be constructed to encapsulate information about individual 
participating biological molecules and information about relationships between them. 
For example, a molecular relational graph encapsulating gene expression data defines 
vertices as genes and edges as connections between genes with significantly correlated 
expression profiles. In another example, a molecular relational graph representing 
metabolic pathway defines vertices as metabolite molecules, edges as connections 
between metabolites related to each other by a single biochemical reaction, and edge 
weights as enzyme that catalyze the reaction between the connected metabolites. As 
used herein, the terms "graph", "graphing", "graphical" are intended to refer to 
mathematical representations recognized as graphs and are not intended to be limited to 
be limited to visual depictions of data (although such visual depictions of data are 
encompassed by the disclosed method). 

Possible types of molecular relational graph include but are not limited to the 
following: 

molecular relational graph representing physical mapping of genes, open 
reading frames, single nucleotide polymorphisms, expressed sequence tags, sequence 
tag sites, or a combination thereof; molecular relational graph representing genetic 
mapping of genes, open reading frames, single nucleotide polymorphisms, expressed 
sequence tags, sequence tag sites, or a combination thereof; molecular relational graph 
representing radiation-hybrid mapping of genes; molecular relational graph 
representing orthologous relationships between genes; molecular relational graph 
representing paralogous relationships between genes; molecular relational graph 
representing homologous relationships between genes; molecular relational graph 
representing structural relationships between proteins; molecular relational graph 
representing gene expression regulation; molecular relational graph representing gene 
translation regulation; molecular relational graph representing protein-protein 
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interactions; molecular relational graph representing protein-DNA interactions; 
molecular relational graph representing enzyme functions; molecular relational graph 
representing chemical metabolism; molecular relational graph representing cellular 
signal transduction pathways; and molecular relational graph representing functional 
5 gene annotation, functional pathways, functional groups, or a combination. 

Illustration 7: Construction of a molecular relational graph representing 
gene expression data. 

Microarray technique has been used widely to measure expression patterns for 
thousands of genes simultaneously. This technique provides a powerful approach for 
10 characterizing gene functions in whole-genome scale. In a typical experiment, 
microarray measurements of gene expression are performed under multiple 
experimental conditions or at multiple time points of a temporal biological process. 
The expression profiles of genes across the treatment are then compared and analyzed. 
The analyses usually consist of a quantification and/or classification of genes into those 
\Q 15 that display similar expression profiles across the experimental conditions. For 

a _ : 

example, if the experimental conditions consist of different time-points in a biological 
process, degree of temporal correlation of expression level for different genes is seen to 
quantify probability of co-regulation of the genes. 

A molecular relational graph representing co-regulation of genes can be 
O 20 constructed by, for example, defining vertices to represent the genes. The method for 

defining a vertex representing a gene is described in Illustration 3. In this type of 
graph, an edge connecting a pair of vertices represents the transcriptional co-regulation 
relationships between a pair of genes represented by the vertex pair. Using methods 
described in Illustrations 4-6, an edge in this type of graph can include following 
25 information items: 

(1) Labels of input and output vertices representing the two genes - 
vertex_labell and vertex_label2. 

(2) Assignment to variable DIRECTED dependent on experiment. 

(3) An edge weight table contains (property _name, value) pairs such as: 
30 (Relationship_type, co-regulation of expression) 

(Relationship_measurement, Pearson's correlation coefficient) 
(Pearson's correlation coefficient, 0.9). 
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As an example, a molecular relational graph representing microarray 
hybridization data for gene expression during the yeast cell cycle (Spellman et al. 
(1998)) was constructed. Pearson's correlation coefficients for the expression profiles 
of a selected set of gene pairs were computed and used as a metric to measure the co- 
regulation relationship and stored in the edge weight table for the edges connecting 
each pair of genes. The resulting molecular relational graph is a completely connected 
graph in which each vertex is connected to every other vertex. A "threshold" graph- 
operation can be performed on the edges of the graph to produce a less densely 
connected graph depicting only the stronger co-regulated relationships. A threshold 
operator zf G,crit) removes vertices or edges from graph G, dependent on the criterion 
set by a conditional statement <crit>. Figure 4 shows an example where a threshold 
operator was applied to the co-regulated yeast molecular relational graph using <crit> 
- if (correlation < 0.6). This operation reveals the co-regulation of expression 
relationships between genes, graded by a degree of confidence. The degree of 
confidence is determined by the threshold parameter. 

Illustration 8: Construction of a molecular relational graph representing 
gene function data. 

A large amount of knowledge about the functions of genes has been 
accumulated in research and documented in research literature. However, large-scale 
systematic exploration and comparison of this body of knowledge with research data 
such as whole genome gene expression profiling data has been hampered by the lack of 
an annotation system that organizes the knowledge into a form enabling transformation 
of the literature into computable quantities. To overcome this obstacle, Gene Ontology 
is the first of such knowledge representation that transforms a large body of knowledge 
about gene functions into a computable collection of annotations (The Gene Ontology 
Consortium (2000)). In Gene Ontology (GO), a comprehensive set of descriptions of 
gene functions is included in the system and each of these descriptions is assigned a 
unique GO identification number (ID). The descriptions are organized in a way such 
that descriptions of related functions are connected to each other in a hierarchical tree 
structure. This tree structure presents the relations between functional descriptions. A 
gene with known function(s) can be assigned one or more GO IDs. Given functional 
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annotations of genes by GO IDs, the disclosed graphs can be used as an effective 
approach to reveal functional relationships for a large number of genes. 

To create a molecular relational graph based on GO annotations of genes, 
vertices representing all genes of interests can be defined. Vertex definition is 
described elsewhere herein (see, for example, Dlustration 3). An edge in the graph 
connects a pair of vertex and encapsulates functional relationship between the two 
genes represented by the vertex pair. An edge can be defined, for example, by the 
following: 

(1) Labels of input and output vertices representing the two genes - 
vertex_labell and vertex_label2 

(2) Assignment to variable DIRECTED depending on the GO function. 

(3) An edge weight table of properties of the functional relationship stored as 
(property jname, value) pairs. As an example, protein product of gene 2 is a 
transcriptional factor that activates the transcription of gene 1. To represent this type of 
functional relationship, an edge weight table can be constructed to contain 
(property _name, value) pairs such as: 

(Relationship_type, transcriptional regulation) 

(Relationship_measurement, K) 

(K, <transcriptional_activation_rate_constant>). 

K is a rate constant used to characterize the kinetics of transcriptional activation 
process. 

When multiple functional relationships happen between a pair of genes, a graph 
can be constructed for each functional type and merged with the AND graph operator 
as described elsewhere herein. Figure 3 shows an example of using Gene Ontology 
(GO) functional annotations for a selected set of yeast genes. Yeast GO functional 
annotation data were imported from the Web site of Gene Ontology Consortium 
(http://www.geneontology.org/) and used to define edges between the subset of genes. 
Connected genes share the same unique GO functional identifier. The graph in Figure 
3 clearly shows known functional relationships for a subset of yeast genes. More 
importantly, from an inspection of the molecular relational graph, one can deduce 
higher-order functional gene relationships not previously characterized. 
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C. Operators 

Operators used in the disclosed method (referred to herein as operators, 
molecular relational graphing operators, or gene-graph operators) are any operation or 
function that can be used to manipulate, transform, combine, split, separate, filter, or 
otherwise alter one or more graphs to produce one or more product graphs. Operators 
that can be used on the disclosed graphs can manipulate the graphs as objects, much as 
mathematical operators manipulate numbers. Like mathematical operators, molecular 
relational graphing operators and gene-graph operators allow direct manipulation of 
graphs using graph operations such as difference, addition, and intersection. Operators 
can be recursive. The disclosed method is not limited to the operators described herein. 
Numerous graph operators and graph manipulation procedures are known and can be 
used in the disclosed method. As used herein, "operation" refers to the use of one or 
more operators on one or more graphs. The disclosed graphs are generally 
mathematical constructs describing biological molecules that can be manipulated, 
transformed, combined, split, filtered or otherwise altered using any relevant 
mathematical operator. 

Operators are defined for computing molecular biological information using 
graphs defined above as operand(s). Rules can be defined for construction of 
biologically meaningful computations. Two or more graphs can be manipulated to 
yield a third graph. Such manipulations allow synthesis of disparate biological 
information encapsulated in different molecular relational graphs. 

Graph operators include unary operators, binary operators, and n-nary operators. 
Useful unary operators include, for example: 

"Threshold edges" which deletes all edges below or above a particular range of 
edge weights; 

"Threshold vertices" which deletes all vertices below or above a particular range 
of vertex parameters; 

"Subset" which is inclusive of only certain edges or vertices (if applied to 
vertices, inapplicable edges are also deleted); 

"Split" which divides one graph into two graphs; 

"Convert graph" which converts a graph from one type to another so that graphs 
of different types can be comparable. 
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Useful binary and n-nary operators include: 

"And" which, given n graphs, finds the common subset of vertices and edges 
and outputs the graph containing only the common vertices and edges; 

"Or" which, given n graphs, finds the union of all vertices and edges and 
outputs the graph containing the union; 

"Addition" which grafts two different graphs A and B together if the two 
different graphs have common vertices; 

"Subtraction" which deletes from a third graph X any vertices common to a first 
graph A and a second graph B; 

"Filtration" which compares and generates a graph X wherein all edges 
(vertices) in compared graphs A, B, etc. that are not also in X are deleted; 

"Consensus" which provides an X% consensus graph of graphs A, B, etc. which 
is defined as a graph consisting of all vertices and edges present in X% or more of the 
graphs, A, B, etc. 

Useful Vertex and Edge operations used in the present invention include: 
"Delete" which deletes a vertex (edge); 
"Add" which adds a vertex (edge); 

"Combine" which combines two or more vertices into one retaining the edges to 
all other vertices or combines two or more edges into a hyper-edge; 

"Examine vertex" which shows information contained in a vertex such as its 
label (gene name), mapping location, amino-acid composition, and can show, for 
example, information obtained through an outside database via a URL linkage; 

"Examine edge" shows information contained in an edge such as 
activation/repression nature of the gene relationship, catalytic rate constant of the 
enzyme reaction, and binding affinity between two protein molecules. 

Operators can be depicted using symbols. This can aid in combining operators 
into sets and series, and in constructing complex operators. An example of a system of 
operator symbols and their use is described below. Additional operators are also 
provided below. 

1. Unary Operators (A) 

Threshold edges (Ai): Delete all edges below (or above) a particular range of 
edge weights. 
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Threshold vertices (A 2 ): Delete all vertices below (or above) a particular range 
of vertex parameters. 

Subset (A3): Inclusive of only certain edges or vertices. If applied to vertices, 
irrelevant edges are also excluded. 
5 Split (A4): Divide one graph into two graphs. 

Find topological sorting for a set of vertices (A 5 ): Find a linear order for a set of 
vertices in a graph such that any graph traversal path constructed from the sorting 
preserves the original order of vertex-to-vertex connection in the graph. 

Find shortest path from vertex A to B (A*): Identify a path starting from vertex 
10 A and ending at vertex B. The number (if un-weighted graph) or the sum of weights (if 
weighted graph) of edges involved in the path is minimum compared to any other 
possible path. 

Find shortest path between each pair of vertices (A 7 ): Identify a path for each 
pair of vertices. The path connects two vertices in the pair and the number (if un- 
15 weighted graph) or the sum of weights (if weighted graph) of edges involved in the 
CO path is minimum compared to any other possible path. 

Find transitive closure (As): Construct for a graph a vertex reachability matrix 
in which the value of an element located at i-th row and 7-th column represents vertex j 
1=* is reachable from vertex i if the value equals to 1 or else 0. 

20 Find articulation points (A 9 ): Traverse the graph and identify all vertices the 

deletion of which splits the graph into two or more substructures. An articulation point 
usually represents a junction linking multiple pathways or subsystems, for example, a 
gene that participates in multiple biological processes. 

Find strongly connected components (A10): Traverse the graph and identify all 
25 subsets of vertices whose connections to vertices within the same subset are much 

denser than are connections to vertices outside the subset. A subset usually reflects a 
relatively complete and independent functional group of genes participating in a single 
biological process. 

Find minimum-weight spanning tree (An): Construct a tree from a graph so that 
30 the tree contains all the vertices in the graph and the sum of weights of all edges in the 
tree is minimum. A tree is a graph with properties: a) any two vertices are connected 
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by precisely one path; b) no vertex can reach itself through a path including zero or 
more edges and/or vertices. 

Find maximum-weight spanning tree (Ai 2 ): Construct a tree from a graph so 
that the tree contains all the vertices in the graph and the sum of weights of all edges in 
the tree is maximal. 

Find fundamental circuits (A13): Find a set of circuits in a graph so that any 
circuit present in the graph can be derived from a ring-sum of a combination of 
elements in the set. A ring-sum of two graphs d = (Vi, E x ) and G 2 = (V 2 , E 2 ) is the 
graph ((ViuV 2 ), ((E1UE2) - (EinE 2 )). 

Find fundamental cut-sets (Am): Find a set of cut-sets in a graph so that any cut- 
set of the graph can be derived from a ring-sum of a combination of elements in the set. 
A cut-set of a connected graph or component is a set of edges whose removal will 
disconnects the graph or colmponent. 

Find the capacity of a cut-set (A15): Calculate the flow capacity of a cut-set of a 
graph. Given a vertex, x, as the source and another vertex, y, as the sink of a network 
N, a flow for N associates a non-negative integer fiu.v) with each edge (w,v) of N, such 
that for all vertices v, other than x or y: ^ /(«, v) = ^ f(v, u). An edge capacity c(w,v) 

u u 

is defined as the maximum of f(u,v) for the corresponding edge. A cut-set of a graph 
(V, E) partitions vertices into two sets (P, P ) such that Pn P = 0 and Pu P = V. The 
capacity of the cut-set is then defined as ^c(w, v). 

ueP 

Condense graph (Ai 6 ): Collapse each component in a graph into a hyper-vertex 
and replace edges incident to and from the component with edges incident to and from 
the hyper- vertex. 

Convert graph (An): Transform a graph from one type to another so that graphs 
from different sources can be compared. 

Find connected components (Ai 8 ): Identify all connected components in a 

graph. 



31 



2. Binary and n-nary Operators (S) 

AND (Hi): Given n graphs, find the common subset of vertices and edges. 
Output the graph containing only the common vertices and edges. 

OR (S2): Given n graphs, find all vertices and edges. Output the graph 
containing all vertices and edges present in either graph. 

Addition (S3): If two different graphs have common vertices, merge the two 

graphs. 

Subtraction (H4): Given graph A and graph B with common vertices, subtraction 
of graph B from graph A is the operation that deletes from graph A all vertices common 
to graph B, thus producing graph C, such that C= A-B. 

Filtration (S5): A filtration of graphs by some graph X is the process of deleting 
all edges (or vertices) in each graph that are not also present in graph X. 

Consensus (He): An X% consensus graph is the graph consisting of all vertices 
and edges present in X% or more of the graphs on which the operation is performed. 

Isomorphism (H7): Given graphs Gi = (Vi, Ei) and G2 = (V2, E2), find a graph 
G 3 = (V 3 , E 3 ) such that: a) there is a bijection/i: Vi S — » V 3 such that {/i(x)/ 2 (y)}eE 3 if 
and only if {x,y}eEi; b) there is a bijection/2* V 2 S — > V 3 such that {/20O/2(y)}eE 3 if 
and only if {x,y}eE 2 where Vi S and V2 S are subsets of Vi and V 2 respectively. A 
bijection is a function f: A — > B if it is both an injection (one-to-one) and a surjection 
(the reverse is also one-to-one)(Ore, Theory of Graphs, American Mathematical 
Society, Providence, RI (1962)). 

3. Vertex and Edge Operators OF) 
Delete C¥\): Remove a vertex (or edge). 
Add OF2): Insert a vertex (or edge). 

Union OF3): Combine two or more vertices into one vertex retaining the 
previously existing edges to all other vertices. Combine two or more edges into a 
hyper-edge. 

Disassemble OF4): Disassemble a hyper-vertex and/or a hyper-edge formed as a 

result of Union operation into original set of vertices and/or edges. 

Examine vertex OF5): Show information contained in a vertex, such as its label, 

gene name, mapping location, amino-acid composition, and URL to external databases. 
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Examine edge C¥e): Show information contained in an edge such as 
activation/repression nature of the gene relationship, catalytic rate constant of the 
enzyme reaction, or binding affinity between two protein molecules. 

4. Rules 

Any computation on molecular relational graphs using molecular relational 
graph operators can be constructed by following rules. The following rules are 
examples of useful rules. In the rule definitions, Gi, G 2 , G3, . G n and G each 
represents a different molecular relational graph and 0 is an empty set. 

(i) Rules of modifiers 

Rules of modifiers can define the syntax for using modifier-style operators, A 
and W. An operator of this type operates on a single input graph: 
Ai(G0 = G 2 , where i = { 1, 2, 3, 6, 7, 11, 12, 16, 17} 
Ai(G) = S, where S = {Gi, G 2 , ...} andi = {4, 10, 13, 14, 18} 
Ai(G) = S, where S = { Vi, V 2 , ... } and i = {5, 9} 
Aj(G) = M, where M is a reachability matrix and i = {8} 
Ai(G) = C, where C 6 Rand i = { 15 } 
Y(G, S) = G, where S = {VI, V2, ... } and i = {8} 

(ii) Rules of binary operation 

Rules of binary operation can define the syntax for using binary operators, 
which take two input graphs and produce an output graph: 
Gi Hi G 2 = G 3 , where i - { 1, 2, 3, 4, 5, 7} 

(iii) Rules of n-nary operation 

Rules of n-nary operation can define the syntax for using n-nary operators, 
which take more than two graphs as input and produce different types of output: 
Hi (Gi, G 2 , G 3 , . . G n ) = G, where i = { 1, 2} 

Hi (S, G) = S\ where S = {Gi, G 2 , G 3 G n }, S* = { Gi\ G 2 \ G 3 ', G n '} and 

= {5}' 

Hi (X%, Gi, G 2) G 3 , G n ) = G, where X% G Rand i = {6} 

(iv) Empty graph laws 

Empty graph laws can define the result of computation for various operators 
when an empty set, 0, is involved in the input: 
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Ai(0) = 0, where i = {1,2,3,4,5,6, 7, 9, 10, 11, 12, 13, 14, 16, 17} 

Ai(0) = M, where M is a reachability matrix with all elements equals to 0 and i 

= {8} 

Ai(0) = O, where i = {15} 
W(0, S) = 0 
^(G, 0) = G 

G Hi 0 = 0, where i = { 1, 6, 7} 

G Si 0 = G, where i = {2, 3, 4, 5} 

Si(0,G 2 ,G 3 , ...,G„) = 0, where i = {1} 

Si (0, G 2 , G 3 , G n ) = Si (G 2 , G 3 , G n ), where i = {2} 

Si (S, 0) = S, where S = {Gi,G 2 , G 3 , ...,G„} andi = {5} 

Si (C, 0, G 2 , G 3 , . . ., G„) = 0, where C G Rand i = {6} 

(v) Idempotency laws 

Idempotency laws can define the result of computation for binary and n-nary 
operators when identical graphs are taken as the input: 
G Si G = G, where i = { 1, 2, 3, 7} 
G Si G = 0, where i = {4, 5} 
3i(G, G, G, ...,G) = G, where i = {1,2, 3,7} 
Si (G, G, G, G) = 0, where i = {5} 

(vi) Commutative laws 

Communitive laws state that, in consecutive binary operations, operands 
involved can exchange positions freely without affecting the end result: 
Gi Hi G 2 = G 2 Hi Gi, where i = { 1, 2, 3, 4, 7} 

(vii) Associative laws 

Associative laws state that the order of a sequence of operations performed by 
binary or n-nary operators can be rearranged without affecting the end result: 
(Gi Hi G 2 ) Hi G 3 = Gi Hi (G 2 Hi G 3 ), where i = { 1, 2, 3, 4, 5, 6, 7 } 

(viii) Distributive laws 

Distributive laws state that the product of a first binary or n-nary operation on 
the product of a second binary or n-nary operation on some objects will yield the same 
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result as the second binary or n-nary operation on the products of the first binary or n- 
nary operation on each of the objects: 

Gi Si (G 2 Sj G 3 ) = (G] Hi G 2 ) Ej (Gi Si G3), where i = { 1, 4, 5, 6, 7}, j = { 1, 2, 

3,4, 6,7},andi*j 

Ai (Gi Sj G 2 ) = (Ai(GO) Sj (Ai(G 2 )), where i = {1, 2, 3, 6, 7, 11, 12, 16, 17}, j = 
{1,2,3,4, 6,7} 

5. Methods for assimilating disparate molecular biological data 

(i) Integration of disparate data sets 

Two or more non-overlapping data sets, {Gi, G 2 , G3, . . ., G n }, can be 
synthesized into a single data set, G: 

G = H 2 (Gi, G 2 , G 3 , . . . , G n ) or G = Gi S 2 G 2 

Two or more overlapping data sets, {Gi, G 2 , G 3 , . . G n }, can be synthesized 
into a single one, G: 

G = S 3 (Gi, G 2 , G 3 , . . G n ) or G = Gi S 3 G 2 

(ii) Filtration of a data set using another data set 

Subtraction of data found in one data set, G 2 , from another data set Gi and yield 
a third data set, Gi f : 

Gi' = Gi H4 G 2 

Filtering out consensus data between one data set, GI, and another data set, G2, 
from data set GI and yield a third data set, Gi': 
d'^Gi H 5 (Gi Si G 2 ) 

(iii) Identification of consensus data from disparate data sets 
Identification of consensus data, G, between two data sets, Gi and G 2 , without 

having to preserve the relationships between biological molecules in original data sets: 
G = Gi Si G 2 

Identification of consensus data, G, between two data sets, Gi and G 2 , such that 
the original relationships between biological molecules are preserved in the resulting 
data set: 

G = Gi S6 G 2 

Identification of consensus data, G, among many data sets, Gi, G 2 , G 3 , G n , 
such that the consensus data appears in more that X% of total number of data sets: 
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G = S 6 (X%,Gi,G 2 ,G 3 , ...,G n ) 

(iv) Identification of unique data for individual disparate data sets 

Identification of data, (Gi, un i q ue, G 2j unique, G 3 , unique, G n , unique,), unique for 
individual data sets, (Gi, G 2 , G 3 , G n ) - method (I): 

5 Gconsensus = Si(Gi, G 2 , G3, . . ., G n ) 

Gi t unique = Gi —4 G CO nsensus 
G 2 , unique ~ G 2 £14 G C onsensus 
G 3) unique = G 3 C*4 G CO nsensus 

* * * 

10 G - G ~aCt 

±KJ ^n, unique ~~ fc — «4 ^consensus 

Identification of data, (Gi, un iq Ue , G 2 , un i que , G 3j unique, • . . , G n , unique,), unique for 
individual data sets, (Gi, G 2 , G 3 , G n ) - method (II): 



Gconsensus = (... ((Gi H 7 G 2 ) E 7 G 3 ) E 7 . . .) H 7 G n 

^ Gi f unique = Gi ^4 Gconsensus 

yj 15 G 2) unique = G 2 £4 Gconsensus 

G 3j unique = G 3 £-4 Gconsensus 



3-J 

lAff Gn, unique = Gn »Sl4 Gconsensus 

I ***** 



(v) Identification of common biological pathways revealed by two 
^ 20 different data sets 

To find a set of biological pathways, S, that are revealed in both data sets, GI 
and G2, one identifies strongly connected components in both graphs first. Then 
condenses those components into hyper- vertices. An isomorphic sub-graph, G, of Gj 
and G 2 is subsequently identified. Pathways can then be isolated from G and stored in 
25 S: 

G = (A 16 (G lf Aio(Gi))) H 7 (A 16 (G 2 , A 10 (G 2 ))) 

S = AisCG), where S is a set of graphs, each of which represents a pathway 
common to both data set Gi and G 2 
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(vi) Identification of biological molecules critical for multiple 

biological pathways 

To identify biological molecules critical for multiple biological pathways (Gi, 
G 2 , G 3 , . . G n ), one identifies articulation points in each graphs first (Vi, V 2 , V 3 , 
V n ) and subsequently find an intersection set, V, of vertex set (Vi, V 2 , V3, . . ., V n ): 

V 1 =A 9 (G 1 ) 

V 2 = A 9 (G 2 ) 

V 3 = A 9 (G 3 ) 

• • * 

V n = A 9 (G„) 

V = V!nV 2 nV 3 n...nV n 

6. Ancillary Functions 

"Find articulation points" which traverses the graph and identifies all the 
vertices that, when deleted, can split graph into two or more substructures; an 
articulation point usually represents the cross-linking point among multiple pathways or 
subsystems, for example, a gene functions in multiple biological processes. 

"Find strongly connected components" which traverses the graph and identifies 
all subsets of vertices whose connections to vertices within the same subset is much 
denser than to the outside vertices; a subset usually reflects a relatively complete and 
independent functional group of genes participating in a single biological process. 

7. Assimilating disparate molecular biological data 

Large-scale and high throughput biological experiments such as whole genome 

gene expression and protein translation profiling produce disparate data of large size. 

The complexity of the relationship information embedded in these data made analysis 

difficult using prior methods. Moreover, these data contain different types of 

relationship information depending on the design and the purpose of the experiments 

generating the data. The heterogeneity of these data presented a serious challenge to 

the integration of information using prior methods. The disclosed method is 

particularly apt for handling the complexity and heterogeneity of data and is thus 

capable of facilitating the integration and understanding of large-size heterogeneous 

biological data. Two examples of the application of the disclosed method to complex 

data are described below and illustrate these capabilities. 
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Illustration 9: Integration of gene expression data with Gene Ontology data 

Microarray gene expression data contain information about expression profiles 
for a large number of genes. From this type of data, gene functions can be inferred by 
comparing expression profiles between genes. Genes having similar expression 
profiles are considered to have high probability of being co-regulated by the same 
transcriptional control mechanism and thus may contribute to the creation of the same 
phenotype. While analyses of newly generated data using state-of-the-art technology 
give tremendous insights into gene functions, discoveries made in previous research 
also accumulate a large body of knowledge that needs to be merged together with 
current progress in order to facilitate the formation of a comprehensive understanding 
of gene functions. One good example of such previously accumulated knowledge is 
Gene Ontology annotations. Integration of gene co-regulation information with 
functional annotation of genes is needed to produce a comparison of these two bodies 
of information. This integration can be done by the synthesis of information 
represented by the disclosed methods. Gene expression data (Spellman et al. (1998)) 
and GO annotation for yeast genes were chosen to illustrate the ability of graph- 
operators to derived integrated representation of heterogeneous information. 

A graph of gene expression profiles was generated from the data as described in 
Illustration 7. In this graph, relationships of expression co-regulation between genes 
are captured by the edges. A second molecular relational graph representing GO 
annotation of genes is generated as described in Illustration 8. To simplify the 
computation, the graph representing GO functional relationships was created as an 
unweighted graph by omitting the step of creating an edge weight table. Since the 
graph of GO functional relationships was an unweighted graph, while the graph of gene 
expression was a weighted graph in which the edge weights were the correlation 
coefficients, the unary operator "convert", c(G, ti, t 2 ), was used to transform a graph 
(G) from one type (ti) to another (t 2 ), so that graphs from different sources can be 
compared. Thus the operator c(G, ti, t 2 ), where ti = WEIGHTED and t 2 = 
UNWEIGHTED, transformed the weighted graph shown in Figure 4 to an unweighted 
graph. 

To integrate the two types of information, the graph of the complete set of GO 
functional relationships (not shown) and a graph of gene expression data (Figure 4) 
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were input to the graph operator "AND". The binary operator "AND" synthesizes 
information from two or more graphs by finding the subset of common edges and 
vertices. The resulting consensus information is shown in Figure 5A. Because only a 
subset of the 6,000+ yeast genes is used to generate Figure 4, the results shown in 
Figure 5A are for illustrative purposes only, and do not represent an exhaustive survey. 
Figure 5A shows two connected component structures representing two distinct sets of 
genes. These sets represent those genes whose GO functional relationships are 
concordant with their expression pattern relationships. 

Illustration 10: Exploratory thresholding of gene expression data. 

In a weighted graph representing co-expression relationships of genes, every 
vertex can be connected with all other vertices through edges. The edge-weights, 
correlation coefficients, for this type of graph quantifies the degree of co-expression. 
The quantitative information in the correlation coefficients can be used to generate a 
coarser representation graph showing only those relations with high confidence. For 
this purpose, the edge filtering operation on molecular relational graphs can be 
performed by the "threshold" operator r(G,crit), which removes vertices or edges from 
graph G, dependent on the criterion set by a conditional statement <crit>. 

As an example of exploratory thresholding applied to gene expression graphs, 

threshold operations were performed on the graph shown in Figure 4 to determine 

whether stronger correlations in gene expression are related to functional relationships. 

That is, it was asked whether the structure shown in Figure 5A can be recovered from 

the graph shown in Figure 4 alone by including only the strongest co-expression 

relationships. In fact, both of the connected graph components seen in Figure 5A 

appear in gene expression graphs thresholded at 0.9 (Figure 5B), 0.8 (Figure 5C), and 

0.7 (Figure 5D). Higher-stringency thresholding produces fewer gene-relationship 

structures in the expression data, but more of the structures produced are supported by 

the GO functional annotations. This suggests a quantitative relationship between 

concordant expression of genes and their functional interaction. In addition, Figure 5 

shows that the expression data also imply some gene relationships (marked by V in 

Figure 5B, 5C, and 5D) which are not apparent in the GO annotation graph (Figure 3). 

Careful examination shows that a higher-order relationship documented in the GO tree 

can account for these expression relationships (Figure 5E). This exercise demonstrates 
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how a novel functional inference could be made through the power of integrative 
analysis using the disclosed method. Operations used to generate Figure 5 are 
summarized in the Table 4. 

Table 4. Operations used to generate the molecular relational graphs shown in 

Figure 5. 



Graph A 


Graph B 


Operator 


Resulting 
Graph 


GO graph 


Gene expression graph 


AND 


Figure 5A 




Gene expression graph 


T(G,crit) 

<crit> = if (correlation < 0.9) 


Figure 5B 




Gene expression graph 


l(G,crit) 

<crit> - if (correlation < 0.8) 


Figure 5C 




Gene expression graph 


T(G,crit) 

<crit> = if (correlation < 0.7) 


Figure 5D 



D. Implementation 

In one embodiment, a software program for GGO can be developed using the 
JAVA programming language. This program has two principal features, the first being 
the implementation of molecular relational graph objects and the ability to persist to a 
local database, and the second being implementation of the set of operators that can be 
performed on the gene-graphs. This software performs the task of integrating the data 
from microarray gene expression analysis, Gene Ontology annotation, and protein- 
protein interaction analysis into a GGO data model functionalities for pathway analysis, 
critical gene identification, gene-action subsystem identification, and pathway 
comparison. Since the molecular relational graphing model is best illustrated using a 
graphical approach, in a preferred embodiment, the software provides visualization 
essential for the demonstration of the data resulting from the computation using GGO 
data model. In a preferred embodiment, the visualization software is based on three 
development resources: JAVA 2D and JAVA3D API libraries developed by SUN 
MICROSYSTEM which provide classes for writing two- and three-dimensional 
graphics applications; Open source software Graphviz developed by AT & T 
Laboratory (www.research.att.com/sw/tools/graphviz/) which is a set of tools for 
construction and geometric presentation of graphs and networks with a publicly 
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available source code allowing use to build complex visualization functionality; and 
commercially available graphics API libraries developed by Advanced Visual Systems. 

Standard analysis techniques can be integrated into this analysis platform by 
incorporating standard commercial software packages. This allows the system to use 
many analysis features, such as clustering analysis, from other packages for preliminary 
data processing. The resulting data is then ported into the molecular relational graphing 
model for high-level analysis. 

An Unified Modeling Language entity diagram of GGO objects employed in the 
design of this software is depicted in Figure 15. 

The analysis capability of the molecular relational graphing data model is 
exemplified in part by the following conversion of genomic information into graph 
structure. Software has been developed to convert genomic information to graph 
structure. Various graph operators have also been implemented for the MRG model, 
including, but not limited to, add and delete vertex, add and delete edge, threshold 
edges, subset, graph AND, and graph OR. Using these programs, data from microarray 
gene expression assays, protein-protein interaction assays, and Gene Ontology 
functional annotation have been encoded into graph structures. Further, a set of graph 
visualization tools have been incorporated into the program. 

Exemplary results are shown in Figures 2 through 5. In Figure 2, data were 
imported from the analysis of the yeast (Saccharomyces cerevisiae) genome and 
encoded into gene-graphs. In this application, 1,004 genes and 957 protein-protein 
interactions documented in Uetz et al. (2000) were graphed. The resulting visualization 
reveals structural complexities such as the subset of strongly connected components 
seen in the middle of Figure 2. 

Similarly, Figure 3 shows a graphical representation of functional relationships 
found in the Gene Ontology (GO) database for a selected set of yeast genes. The 
resulting graph encapsulates previous knowledge of the function of these genes. A 
comprehensive view of the functional relationships among the genes is clearly revealed 
by the gene-graph. Importantly, the gene-graph representation reveals higher-order 
functional gene relationships not previously characterized. 

Quantitative relational data such as correlations can also be represented as a 
graph structure. As an example of this, microarray hybridization data were analyzed 
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for gene expression during the yeast cell cycle (Spellman et al. (1998)). The expression 
profile correlations of all gene pairs were computed and used as a metric to define the 
edge weight for the edges connecting each pair of vertices, here defined as genes. The 
gene-graph thus generated encapsulates the relationships of the gene expression 
profiles. The unary operation "thresholding" converts quantitative relational 
information into more intuitive qualitative information with a tunable parameter. A 
threshold operation on the graph of gene expression was performed. A threshold of 0.4 
was chosen, where a value of 0 corresponds to no correlation, and a value of 1 to 
complete correlation. In this threshold operation, edges were deleted if their weights 
were greater than or equal to 0.4. The resulting graph is shown in Figure 4. This 
operation reveals the expression relationship between genes, graded by the degree of 
confidence as measured by a quantitative parameter. 

Information from two or more kinds of gene-graph can be synthesized using the 
graph operation AND. Figure 5 presents such a synthesis of information between the 
functional relationship indicated by the GO gene-graph and the Spellman et al. 
expression study. The AND operator was used with different threshold operators on 
the expression graph to demonstrate how graph operators can be combined to yield a 
flexible set of information syntheses. Figure 5 A, shows the results of an AND 
operation between the GO annotation graph and gene expression graph thresholded at 
the 0.4 level. The result produces two connected component structures representing 
two distinct sets of genes whose functional relationships are concordant with their 
expression pattern relationships. Both structures appear in expression gene-graphs 
thresholded at 0.1 (Figure 5B), 0.2 (Figure 5C), and 0.3 (Figure 5D). Higher- 
stringency thresholding produces fewer gene-relationship structures in the expression 
data, but more of the produced structures are in conformity with the GO data. This 
indicates a quantitative relationship between concordant expression of genes and their 
functional interaction. Figure 5 shows a relationship between genes implied by the 
expression data that is not apparent in the GO data (marked by V). However, careful 
examination shows that a second order interaction documented in the GO accounts for 
the expression relationship (Figure 5E). This is a novel discovery mediated by the 
power of integrative analysis from the GGO model of the present invention. 
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Accordingly, as demonstrated herein, gene-graph analysis provides a powerful 
tool for the analysis of large genomic data sets and the discovery of novel gene 
relationships, as well as for the corroboration of relational data by drawing consensus 
from disparate sources of information. Further enrichment of the algorithmic 
5 operations on the gene-graph by adding new theoretical and heuristic components can 
greatly expand the potential of this analytical technique and transform it into a 
significant discovery tool for genome-scale data analysis. 

The disclosed method can be produced and used at varying levels from software 
components to integrated packages with user-interface which allows a wide range of 
10 application. Different graph manipulation tools can be implemented, for example, as 

reusable JAVA components. In addition, GGO software may be readily interfaced with 
other software packages, such as common statistical packages. A useful component of 
the integrative data analysis package of the disclosed method is to enable preliminary 
data processing, such as cluster analysis. Common statistical packages could be used to 
15 provide such analyses. Thus, all or part of the disclosed method can be implemented as 

S"-_ a 

macros and routines to interface statistical analysis packages such as SAS, SPSS, 

m 

^ SPLUS using the GGO data model. 

Software design process for implementing the disclosed method preferably can 
employ the object-oriented notation, UML (Unified Modeling Language, Booch et al.), 
Q 20 to document requirements, classes, class behavior, and class dependencies of molecular 

relational graphing software. A UML entity diagram of a selection of molecular 
relational graphing objects is shown in Figure 15. In order to capture the architectural 
design of the molecular relational graphing software, user interface story-boards, use 
case diagrams, sequence diagrams, and class hierarchy diagrams can be developed. 

25 E. Embodiments 

The disclosed method, structures, and compositions can be further understood 
with the following descriptions of some of their forms and embodiments. 

One embodiment of the disclosed method is a computer-implemented method 
for performing an operation upon one or more graphs, wherein each graph can 
30 represent a set of relationships between a set of biological molecules, wherein each 
graph can comprise vertices representing the biological molecules and edges 
representing the relationships between the biological molecules, where the method 
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comprises performing one or more operations on the one or more graphs to produce one 
or more product graphs. 

Another embodiment of the disclosed method is a computer-implemented 
method for performing an operation upon a graph, where the graph can represent 
relationships between biological molecules and can have vertices representing the 
molecules and edges representing the relationships, where the method comprises 
identifying a subset of zero or more of the edges, identifying a subset of zero or more of 
the vertices, and performing a unary operation upon the identified subset of edges and 
vertices to produce a product graph. As used herein, "identifying a subset" of vertices 
and/or edges refers to selecting, using any desired criteria, those vertices and/or edges 
in a set of vertices, set of edges, and/or graph(s) having or lacking one or more of the 
desired criteria features. 

Another embodiment of the disclosed method is a computer-implemented 
method for representing relationships between biological molecules using one or more 
graphs each having vertices and edges, where the method comprises representing a set 
of biological molecules, wherein each molecule can be represented by a vertex of the 
graph, and representing a set of relationships between the biological molecules, 
wherein each relationship can be represented by an edge of the graph, wherein the edge 
connects two vertices, wherein the graph can be produced by performing one or more 
operations on one or more input graphs to produce the one or more graphs. The 
disclosed graphs represent relationships between biological molecules. 

One embodiment of the disclosed composition is a computer program product 
for performing an operation upon one or more graphs, wherein each graph can 
represent a set of relationships between a set of biological molecules, wherein each 
graph can comprise vertices representing the biological molecules and edges 
representing the relationships between the biological molecules, where the computer 
program product comprises a computer data medium on which is carried a means for 
performing one or more operations on the one or more graphs to produce one or more 
product graphs. 

Another embodiment of the disclosed composition is a computer program 
product for performing an operation upon a graph, where the graph can represent 
relationships between biological molecules and can have vertices representing the 
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molecules and edges representing the relationships, where the computer program 
product comprises a computer data medium on which is carried a means for identifying 
a subset of zero or more of the edges, a means for identifying a subset of zero or more 
of the vertices, and a means for performing a unary operation upon the identified subset 
of edges and vertices to produce a product graph. 

Another embodiment of the disclosed composition is a computer program 
product for representing relationships between biological molecules using a graph 
having vertices and edges, where the computer program product comprises a computer 
data medium on which is carried a means for representing a set of biological molecules, 
wherein each molecule can be represented by a vertex of the graph, and a means for 
representing a set of relationships between the biological molecules, wherein each 
relationship can be represented by an edge of the graph, wherein the edge connects two 
vertices. 

Another embodiment of the disclosed method is a computer-implemented 
method for representing relationships between biological molecules using a graph 
having vertices and edges, where the method comprises representing a set of biological 
molecules, wherein each molecule can be represented by a vertex of the graph, and 
representing a set of relationships between the biological molecules, wherein each 
relationship can be represented by an edge of the graph, wherein the edge connects two 
vertices. 

Another embodiment of the disclosed composition is a representation of 
relationships between biological molecules comprising one or more graphs each having 
vertices and edges, each graph comprising a set of biological molecules, wherein each 
molecule can be represented by a vertex of the graph, and a set of relationships between 
the biological molecules, wherein each relationship can be represented by an edge of 
the graph, wherein the edge connects two vertices, wherein the graph can be produced 
by performing one or more operations on one or more input graphs to produce the one 
or more graphs. 

Another embodiment of the disclosed composition is a data structure 
comprising a representation of relationships between biological molecules, where the 
representation can comprise a graph having vertices and edges, where the graph 
comprises a set of biological molecules, wherein each molecule can be represented by a 

45 



vertex of the graph, and a set of relationships between the biological molecules, 
wherein each relationship can be represented by an edge of the graph, wherein the edge 
connects two vertices. A data structure is any form of data, information, and/or objects 
collected, organized, stored, and/or embodied in a composition or medium. A 
5 molecular relational graph stored in electronic form, such as in RAM or on a storage 
disk, is a type of data structure. 

Another embodiment of the disclosed method is a computer-implemented 
method for graphically representing relationships between biological molecules using a 
graph having vertices and edges, where the method comprises displaying a 
10 representation of a set of biological molecules, where each molecule can be graphically 
represented by a vertex of the graph; and displaying a representation of a set of 
relationships between the molecules, where each relationship can be graphically 
represented by an edge of the graph, where each edge can have an associated 
84 description, wherein the edge connects two vertices. As used herein, a graphical 

15 representation is a visual representation of a graph. 
Jijf Another embodiment of the disclosed method is a computer-implemented 

method for performing an operation upon a graph, where the graph can represent 
relationships between biological molecules and can have vertices representing the 
molecules and edges representing the relationships, where the method comprises 
20 displaying the graph; identifying a subset of zero or more of the edges; identifying a 

subset of zero or more of the vertices; performing a unary operation upon the identified 
subset of edges and vertices; and displaying a product graph resulting from the unary 
operation. 

Another embodiment of the disclosed method is a computer-implemented 
25 method for performing an operation upon a set of n graphs, where each graph can 

represent relationships between biological molecules and can have vertices representing 
the molecules and edges representing the relationships, where the method comprises 
performing an n-nary operation upon the n graphs; and displaying a product graph 
resulting from the n-nary operation. 
30 Another embodiment of the disclosed composition is a computer program 

product for graphically representing relationships between biological molecules using a 
graph having vertices and edges, where the computer program product comprises a 
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computer data medium on which is carried a means for displaying a representation of a 
set of biological molecules, where each molecule can be graphically represented by a 
vertex of the graph; and a means for displaying a representation of a set of relationships 
between the molecules, where each relationship can be graphically represented by an 
edge of the graph, each edge having an associated description. 

In these or other embodiments disclosed herein, the method or composition can 
have any or a combination of the following features. For example, the operations can 
comprise finding a common subset of vertices and edges in a plurality of graphs; 
merging a plurality of graphs having one or more common vertices or edges; deleting 
vertices and edges present in a first graph that are not present in a second graph; 
combining the edges and vertices of a plurality of graphs; finding a common subset of 
vertices and edges present in a predetermined percent of a plurality of graphs; finding a 
common subset of vertices and edges in a plurality of graphs, and deleting the common 
subset of vertices and edges from each of the graphs to produce a plurality of graphs 
each with a unique set of vertices and edges; deleting all edges beyond a selected range 
of edge weights; dividing one graph into two graphs; using an AND operation to find 
the common subsets of vertices and edges of n graphs; or any combination of these 
and/or other operations. Any of the operations can be a recursive operation. 

The set of biological molecules can comprise more than one type of biological 
molecule or can be all of the same type of biological molecule. The biological 
molecules can be, for example, selected from the group consisting of genes, open 
reading frames, expressed sequence tags, single nucleotide polymorphisms, sequence 
tag sites, nucleic acids, DNA, RNA, mRNA, cDNA, proteins, peptides, enzymes, 
metabolites, carbohydrates, exons, introns, cleavage fragments, restriction fragments, 
amino acid modifications, protein domains, DNA or RNA secondary or tertiary 
structures, nucleic acid motifs, protein motifs, and metal ions. 

The set of relationships can comprise more than one type of relationship or can 
be all of the same type of relationship. The relationships can be, for example, selected 
from the group consisting of physical distances between genes, open reading frames, 
single nucleotide polymorphisms, expressed sequence tags, sequence tag sites, or a 
combination thereof; genetic distances between genes, open reading frames, single 
nucleotide polymorphisms, expressed sequence tags, sequence tag sites, or a 
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combination thereof; protein-protein interactions; protein-nucleic acid interactions; 
gene expression regulation; protein expression regulation; cellular signal transduction 
pathways; sequence similarity between genes or proteins; structural similarity between 
proteins; radiation hybrid mapping distances between genes, open reading frames, 
single nucleotide polymorphisms, expressed sequence tags, sequence tag sites, or a 
combination thereof; and metabolic pathways. 

The edges can have a variety of values and features. For example, at least one 
edge can comprise a direction; at least one edge can comprise a boolean value 
indicating the presence or absence of an association between the biological molecules 
represented by the vertices connected by the edge (where, in some embodiments, the 
association can be co-expression, co-regulation, or presence or use in the same 
pathway); at least two of the vertices can represent different types of biological 
molecules; at least two edges can represent different types of relationships between the 
biological molecules represented by the vertices connected by the edges; at least one 
edge can represent a plurality of different types of relationships between the biological 
molecules represented by the vertices connected by the edge; at least one vertex can 
represent a plurality of different biological molecules; at least one edge can comprise 
an edge weight; a subset of edges can be edges beyond a selected range of edge 
weights; or any combination of these and/or other features. 

Where an edge comprises an edge weight, the edge weight can represent a value 
characterizing the relationship represented by the edge (where, in some embodiments, 
the value can be a numerical value; at least one edge can comprise an edge weight table 
comprising the edge weight (where, in some embodiments, the edge weight table 
further can comprise one or more additional edge weights); at least one edge weight can 
comprise an indication of a state; at least one edge weight can comprise a spatial 
distance (where, in some embodiments, the spatial distance can represent a physical 
distance between the biological molecules represented by the vertices connected by the 
edge); at least one edge weight can comprise a kinetic measurement; at least one edge 
weight can comprise a distance metric representing a logical relationship between the 
biological molecules represented by the vertices connected by the edge; at least one 
edge weight can comprise a statistical metric representing a logical relationship 
between the biological molecules represented by the vertices connected by the edge; at 

48 



least one edge weight can comprise a value of fuzzy set membership representing a 
logical relationship between the biological molecules represented by the vertices 
connected by the edge; at least one edge weight can comprise a conditional probability 
(where, in some embodiments, the conditional probability can be the probability of a 
5 causal relationship between the biological molecules represented by the vertices 
connected by the edge); or any combination of these and/or other features. 

The disclosed method and compositions can also comprise hyper-edges and/or 
hyper-vertices. For example, at least one of the graphs can comprise at least one hyper- 
edge (where, in some embodiments, at least one of the operations can convert at least 
10 one hyper-edge to a non-hyper-edge); at least one of the graphs can comprise at least 
one hyper-vertex (where, in some embodiments, at least one of the operations can 
q convert at least one hyper- vertex to a non-hyper- vertex); at least one of the graphs can 

5 comprise at least one hyper-edge and at least one hyper-vertex (where, in some 

K embodiments, at least one of the operations can convert at least one hyper-edge to a 

15 non-hyper-edge, at least one of the operations can convert at least one hyper- vertex to a 
non-hyper- vertex, and/or at least one of the operations can convert at least one hyper- 
edge to a non-hyper-edge and at least one hyper- vertex to a non-hyper- vertex); at least 
one of the operations can convert at least one edge to a hyper-edge (where, in some 
embodiments, the hyper-edge can be formed by combining two or more edges); at least 
S 20 one of the operations can convert at least one vertex to a hyper- vertex (where, in some 

embodiments, the hyper- vertex can be formed by combining two or more vertices; at 
least one of the operations can convert at least one edge to a hyper-edge and at least one 
vertex to a hyper-vertex (where, in some embodiments, the hyper-edge can be formed 
by combining two or more edges and the hyper-vertex is formed by combining two or 
25 more vertices); or any combination of these and/or other features. 

The product graph produced or present in any embodiment of the disclosed 
method or composition can be a graph that is modified relative to the graph on which 
the operation is performed. 

As indicated above, the disclosed methods can be performed using a suitable 
30 computer or other electronic system. In the illustrated embodiment of the invention, 

the methods can be performed using a suitably programmed general-purpose computer 
system such as that illustrated in Figure 14. Persons skilled in the art to which the 
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invention pertains will readily be capable of programming the computer system or 
otherwise providing it with suitable software to implement the above-described 
methods. 

Although the software can be structured in any suitable manner and written in 
any suitable programming languages, it can be conceptually considered to include a 
GGO subsystem 102, and a data mining service broker 104. This software executes in 
the memory 106 of the computer in the manner in which application software 
conventionally executes in such computers. Although GGO subsystem 102 and data 
mining service broker 104 are conceptually illustrated as residing in memory 106 for 
purposes of clarity, persons of skill in the art will recognize that in actual operation they 
may not reside in memory 106 simultaneously or in their entireties. Such persons will 
further understand that many other software elements that typically execute in such a 
computer system, such as operating system software, network communication software, 
software utilities, and other application programs are not illustrated for purposes of 
clarity. 

In addition to memory 106, the computer system can include other suitable 
hardware that is typically included in a general purpose computer, such as a processor 
108, a network interface 1 10, a fixed-medium disk drive 1 12 such as a hard disk drive, 
a removable-medium disk drive 1 14 such as a floppy disk or optical disk drive, and 
input/output interface logic 116. The software elements described that embody a 
system of the present invention can be provided via a program product, such as a floppy 
disk 1 18 on which such elements are recorded. Alternatively, the can be provided via a 
network 120 from a remote site. The software elements can be transferred to disk drive 
1 12 for long-term storage, from where they are used during operation of the system by 
loading them into memory 106 as needed, under the control of processor 108, in the 
manner well-understood in the art. 

The user can interact with the computer system using a mouse 122, keyboard 
124 and video monitor or other display 126 in the conventional manner. Thus, where it 
is described above that the user makes a selection or otherwise provides input in 
response to a displayed menu or other output, such steps can be implemented by using 
mouse 122 and keyboard 124 to provide input in response to information output on 
display 126. Note that descriptions above of outputting graphs for the user refer in the 
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illustrated embodiment of the invention to displaying them on display 126. Although 
not illustrated for purposes of clarity, the graphs can alternatively be output to a printer 
(not shown) or any other suitable output device or sent to a remote system via network 
120. Likewise, graphs can be received from such a remote system via network 120 or 
input via any other suitable input device, such as disk 118. Furthermore, as described 
below, users of remote systems can use the illustrated system for data mining purposes. 

As illustrated in further detail in Figure 6, GGO subsystem 102 can include a 
graph computation manager 130, a graph visualization engine 132, a graph computation 
engine 134 and a graph database 136. Graph computation manager 130 can interface 
not only with graph database 136 but also with other inside databases 140 and outside 
databases 142. Graph computation manager 130 also interfaces with data mining 
service broker 104. The other inside databases can be databases containing 
representations of genes, open reading frames, expressed sequence tags, single 
nucleotide polymorphisms, sequence tag sites, nucleic acids, DNA, RNA, mRNA, 
cDNA, proteins, peptides, enzymes, metabolites, carbohydrates, exons, introns, 
cleavage fragments, restriction fragments, amino acid modifications, protein domains, 
DNA or RNA secondary or tertiary structures, nucleic acid motifs, protein motifs, and 
metal ions. The other inside databases can also contain information about the sample 
collection and experimental processing of the biological materials as captured by a 
Laboratory Information Management System, LLMS. 

Graph computation manager 130 is a middleware component or element that 
performs data mining, visualizes results of data mining, queries previous data mining 
results, and visualizes result data. Graph computation engine 134 is a toolkit/library 
that provides ways to construct graphs and perform graph computations. Graph 
visualization engine 132 creates graphics objects from graph data objects. 

Data mining service broker 104 is a middleware component that communicates 
with a data mining service client 100, decomposes data mining request objects, 
dispatches requests to appropriate subsystems, and receives computational or database 
querying result objects and sends them to data mining service client. 

As illustrated in Figure 7, data mining service client 100 can include a graphical 
user interface (GUI) 150, a request constructor 152, a result unbundler 154, and a 
communications interface 156. 
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As illustrated in Figure 8, data mining service broker 104 can include a client 
manager 160, a client queue 162, a request dispatcher 164, a result dispatcher 166, and 
communications interfaces 167, 168, and 169. 

As illustrated in Figure 9, graph computation manager 130 can include a job 
manager 170, a job queue 172, a graph computational organizer 174, an outside 
database query engine 176, an other inside database query engine 178, a graph database 
engine 180, a graph visualization unit, and communications interfaces 184, 185, 186, 
187, 188, and 189. 

As illustrated in Figure 10, graph computation engine 134 can include graph 
computation engine 190, which can include graph computation executor 192 and graph 
computation library 194, and communications interface 196. 

As illustrated in Figure 11, graph visualization engine 132 can include a graph 
visualization constructor 200 and a communications interface 202. Tom Sawyer GLT 
3.1, referred to in Figures 6 and 11, is only an example of graphical representation 
software that can be used in the graph visualization engine. 

As illustrated in Figure 12, graph computation library 194 can include gene 
graph operator 196, which can include strict graph 198. 

As illustrated in Figure 13, data interface 210 can include a data receiver 212, a 
data transformation engine 214, a request transformation engine 216, and a data 
dispatcher 218. 

Examples 

An example of the disclosed method involving a molecular relational graph of 
genomics data has been implemented using the Java programming language. Software 
has been developed to convert genomics information to graph structure. Using the 
programs, data from microarray gene expression assays, protein-protein interaction 
assays, and Gene Ontology functional annotation (Gene Ontology consortium, 1998) 
have been encoded into graph structures. A set of graph visualization tools is 
incorporated into the programs. 

Data was imported from the analysis of the yeast (Saccharomyces cerevisiae) 
genome, and these data were encoded into molecular relational graphs. As shown in 
Figure 2, the 1,004 yeast genes and 957 protein-protein interactions documented by 
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Uetz et al. (2000) have been graphed. The resulting graph shows structural 
complexities, such as the subset of strongly connected components seen in the middle 
of Figure 2. Similarly, for another data set, data derived from the Gene Ontology (GO) 
annotation for functional relationships of a selected set of yeast genes was encoded. 
The graph shown in Figure 3 was generated by connecting genes that share the same 
unique GO functional identifier. This graph clearly shows known functional 
relationships of the yeast genes. More importantly, from inspection of the molecular 
relational graph, higher-order functional gene relationships not previously characterized 
can be deduced. 

Quantitative relational data such as correlation coefficients also can be 
represented in graph form. Microarray hybridization data for gene expression during 
the yeast cell cycle (Spellman et al, 1998) was analyzed. The correlation coefficients 
for the expression profile of a selected set of gene pairs were computed and used as a 
metric to define the edge weight for the edges connecting each pair of genes. The 
resulting molecular relational graphing (not shown) is a completely connected graph in 
which each vertex is connected to every other vertex. The edges of this graph are 
weighted by the correlation coefficients. However, a "threshold" operation can be 
performed on the edges of the graph to produce a less densely connected graph 
depicting only the stronger relationships. A threshold of 0.6 was used, where a value of 
0 corresponds to no correlation, and a value of 1 to complete correlation. In this 
threshold operation, edges were deleted if their weights are less than or equal to 0.6. 
The resulting graph is shown in Figure 4. This operation reveals the expression 
relationships between genes, graded by a degree of confidence. The degree of 
confidence is determined by the threshold parameter. 

A strength of the disclosed molecular relational graphing model comes from the 
ability to manipulate and combine graphs. In order to demonstrate this capability, a 
small number of graph operators for the molecular relational graphing data model were 
defined, including add vertex, delete vertex, add edge, delete edge, threshold edges, 
convert graph, subset, graph AND, and graph OR. These operators were implemented 
in the example software. 

The molecular relational graph of the complete set of GO functional 
relationships, and the molecular relational graph of expression data shown in Figure 4 
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were used to illustrate graph manipulations. The graph of GO functional relationships 
is an unweighted graph, while the graph in Figure 4 is a weighted graph, in which the 
edge weights are the correlation coefficients. The unary operator "convert" transforms 
a graph from one type to another, so that graphs from different sources can be 
5 compared. The "convert" operator was used to transform the weighted graph shown in 
Figure 4 to an unweighted graph (not shown). 

The binary operator "AND" synthesizes information from two or more graphs 
by finding the subset of common edges and vertices. The "AND" operator was applied 
to the complete set of GO functional relationships (not shown) and the molecular 
10 relational graph of a subset of data from the expression study of Spellman et al. (1998), 
(shown in Figure 4). Figure 5 A depicts this synthesis of information. Because only a 
subset of the 6,000+ yeast genes was used to generate Figure 4, the results shown in 
Figure 5A are merely illustrative, and do not represent an exhaustive survey. Figure 5 A 
shows two connected component structures representing two distinct sets of genes. 
15 These sets represent those genes whose GO functional relationships are concordant 
35 with their expression pattern relationships. 

Additional threshold operations were used on the graph in Figure 4 to determine 

^ whether stronger correlations in gene expression are related to functional relationships. 

M 

That is, it was asked whether the structure shown in Figure 5A can be recovered from 
20 the graph shown Figure 4 alone by subsetting only the strongest pattern relationships. 
Both of the connected components seen in Figure 5 A appear in expression molecular 
relational graphs thresholded at 0.9 (Figure 5B), 0.8 (Figure 5C), and 0.7 (Figure 5D). 
Higher-stringency thresholding produces fewer gene-relationship structures in the 
expression data, but more of the structures produced are supported by the GO data. This 
25 suggests a quantitative relationship between concordant expression of genes and their 
functional interaction. In addition, Figure 5 shows that the expression data also imply 
some gene relationships (marked by V in Figures 5B, 5C, and 5D) which are not 
apparent in the GO molecular relational graph (Figure 3). Careful examination shows 
that a higher-order relationship documented in the GO tree can account for these 
30 expression relationships (Figure 5E). This exercise demonstrates how a novel inference 
can be made through the power of integrative analysis using the disclosed molecular 
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relational graphing data model. Operations used to generate Figure 5 are summarized 
in Table 4. 

Table 5. Operation used to generate the molecular relational graphs shown in Figure 5. 



Graph A 


Graph B 


Operator 


Resulting Graph 


GO graph 


Expression graph 


AND 


Figure 5A 




Expressiongraph 


Threshold at 0.9 


Figure 5B 




Expression graph 


Threshold at 0.8 


Figure 5C 




Expression graph 


Threshold at 07 


Figure 5D 



In summary, the disclosed molecular relational graphing provides a powerful 
tool for the analysis of large genomic data sets and for the discovery of novel gene 
relationships. In addition, it provides an elegant method for the corroboration of 
relational data by drawing consensus from disparate sources of information. Further 
enrichment of the algorithmic operations on the molecular relational graph by adding 
new theoretical and heuristic operators can greatly expand the potential of this 
analytical technique, and transform it into a significant discovery tool for genome-scale 
data analysis. 
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It is understood that the disclosed invention is not limited to the particular 
methodology, protocols, and reagents described as these may vary. It is also to be 
understood that the terminology used herein is for the purpose of describing particular 
embodiments only, and is not intended to limit the scope of the present invention which 
will be limited only by the appended claims. 

It must be noted that as used herein and in the appended claims, the singular 
forms "a "an", and "the" include plural reference unless the context clearly dictates 
otherwise. Thus, for example, reference to "a host cell" includes a plurality of such 
host cells, reference to "the antibody" is a reference to one or more antibodies and 
equivalents thereof known to those skilled in the art, and so forth. 

Unless defined otherwise, all technical and scientific terms used herein have 
the same meanings as commonly understood by one of skill in the art to which the 
disclosed invention belongs. Although any methods and materials similar or equivalent 
to those described herein can be used in the practice or testing of the present invention, 
the preferred methods, devices, and materials are as described. Publications cited 
herein and the material for which they are cited are specifically incorporated by 
reference. Nothing herein is to be construed as an admission that the invention is not 
entitled to antedate such disclosure by virtue of prior invention. 

Those skilled in the art will recognize, or be able to ascertain using no more 
than routine experimentation, many equivalents to the specific embodiments of the 
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invention described herein. Such equivalents are intended to be encompassed by 
following claims. 
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